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Abstract 



, Threshold feedback policies are well known and provably rate-wise optimal selective feedback tech- 

niques for communication systems requiring partial channel state information (CSI). However, optimal 
selection of thresholds at mobile users to maximize information theoretic data rates subject to feedback 
^ I constraints is an open problem. In this paper, we focus on the optimal threshold selection problem, and 

provide a solution for this problem for finite feedback systems. Rather surprisingly, we show that using 
the same threshold values at all mobile users is not always a rate-wise optimal feedback strategy, even 
' for a system with identical users experiencing statistically the same channel conditions. By utilizing the 

00 

I theory of majorization, we identify an underlying Schur-concave structure in the rate function and obtain 

, sufficient conditions for a homogenous threshold feedback policy to be optimal. Our results hold for 

m 

i most fading channel models, and we illustrate an application of our results to familiar Rayleigh fading 

channels. 



I. Introduction 

Consider a multiuser communication system in a fading environment. The channel changes over time, 
and the goal of the base station (BS) is to maximize downlink data rates by taking channel variations 
into account. The BS has multiple transmission antennas, and therefore can possibly communicate with 
a selected subset of users through multiple traffic flows simultaneously. In this setting, selecting the 
users with the best instantaneous signal-to-interference -plus-noise-ratios (SINR) for communication is a 
simple communication strategy that is heuristically expected to maximize downlink data rates. Indeed, 
this is the classical opportunistic communication approach utilizing multiuser diversity to take advantage 
of changing channel conditions for rate maximization H], lEl- 
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For vector broadcast channels, it is well known that such an opportunistic communication strategy 
is asymptotically optimal ||3l. Furthermore, this approach only requires mobile users to feed back their 
partial CSI (i.e., their SINRs), as opposed to other approaches requiring perfect CSI |@], Q. However, 
regardless of partial or perfect CSI feedback, all users contend for the uplink to communicate their CSI 
to the BS in such a setting. This is an impractical burden on the uplink (for large numbers of users), and 
a significant waste of communication resources, which is created by the feedback from users having no 
realistic chance of being scheduled for communication. These considerations motivate the use of selective 
feedback techniques in which only the users with channels good enough are allowed to feed back Q. As 
shown in our companion paper ||7], threshold feedback policies lISl- lfTTI form a rate-wise optimal subset 
of selective feedback policies for downlink throughput maximization subject to constraints on the average 
number of users feeding back. According to a threshold feedback policy, each user decides to feed back 
or not by comparing a metric that captures channel quality (mostly SINR values) with her threshold 
value. Different users are allowed to have different thresholds if such an heterogeneity in thresholds 
maximizes downlink throughput. This paper studies the important issue of how to set threshold values 
at mobile users optimally to maximize downlink data rates without violating finite feedback constraints. 

To this end, we focus on the general setting of vector broadcast channels operating based on the concept 
of opportunistic beamforming [3 |: Each user calculates the SINR on each beam (random orthonormal 
beams in this context) and feeds back the SINR value and the index of the best beam to the BS. Upon 
retrieval of this information, the BS schedules the user with the best SINR on each beam. Certainly, 
threshold feedback policies provide further feedback reductions for such systems, and therefore the 
following question is of practical and theoretical importance to answer: What is the optimal assignment of 
thresholds to users so that aggregate data rates over multiple beams are maximized subject to constraints 
on the average number of users feeding back per beam? 

Recently, a rather surprising answer for this question was presented in lITOl . For a system with identical 
users experiencing the same channel conditions statistically, it is intuitively expected that using the same 
threshold value at all users is rate-wise optimal. However, Inaltekin et. al. in ifTOll showed that this is not 
always true, and established sufficient conditions on SINR distributions for the optimality of homogenous 
threshold feedback policies. The limitation of |[TOl is that they only focused on two-user communication 
systems. In this paper, we use the theory of majorization |[T2l to extend this result to communication 
systems containing arbitrary number of users. 

The rest of the paper is organized as follows. We formally define threshold feedback policies, and 
formulate the optimal threshold selection problem in Section |II1 Section [III] provides some key concepts 
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from the theory of majorization to be utilized in later parts of the paper. The aggregate rate function 
is analyzed (as a function of threshold values and feedback probabilities) in Section JV] This analysis 
leads to sufficient conditions on the SINR distributions to ensure Schur-concavity of the rate function, 
and thereby the optimality of homogenous threshold feedback policies. We illustrate an application of 
our results to familiar Rayleigh fading channels in Section |Vl and Section |Vl] concludes the paper. 

II. System Model 

Consider a vector broadcast channel with M beams. There are n users with each having a single 
receiver antenna. The SINR on beam m at user i is given by 

I X- |2 

7i,m - ^ ^ 2' 
p + l^k=l,kytm \^i,k\ 

where \Xi^k \ is the received instantaneous signal power on beam k at user i, and ^ is the background noise 
power. The random received signal powers at different users are assumed to be independent and identically 
distributed (i.i.d.). This implies i.i.d. fading processes at different users. Let 7^ = (71,1, 7i,2, ■ • ■ , 7i,M)^ G 
be the SINR vector at user i. Beams are identical. That is, for all i € M, the elements of 7j are 
identically distributed (but dependent) with a common marginal distribution F, where J\f = {1,2, ... , n}. 
We will assume that F is continuous, and has the density / with support M+, which are true for many 
fading models. Users feed back according to a threshold feedback policy. We formally define a threshold 
feedback policy as follows. 

Definition 1: We say T = (7i, 72, • • • , 7^)^ is a threshold feedback policy if, for all i G J\f, there is 
a threshold Tj such that Ti{'Ji) generates a feedback packet containing SINR values {'yi^k}i.^x ^^'^ 
only if 7j^fc > Ti for all k £ Zi CI A4, where Ai = {1,2, ■ ■ ■ , M}. We call it a homogenous threshold 
feedback policy if all users use the same threshold r, i.e., ti = t for all i G M. 

Upon retrieval of feedback information, the BS schedules communication to the user with the highest 
SINR on each beam. We define the truncated SINR on beam m at user i as 7^^^ = 7j,ml{7, „>t,}' ^^d 
let 7j be the truncated SINR vector at user i. This vector contains all useful SINR values to be fed 
back by user i on a positive feedback decision. Then, the ergodic sum rate achieved for a given vector 
of thresholds r is given by 

i?(r) = E[r({7jr=i)] = E 

where r ({7i}"=i) is the instantaneous downlink communication rate. Note that zero rate is achieved on a 
beam if no user requests it. Conditioned on an event A (or, a random variable), we define the conditional 



M 



J] log 



1 + max 7j m 

l<j<n 
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rate as R {t\A) = E [r ({7j"=i) |^]. We also let i?™ (r) and ({7j"=i) be the ergodic sum rate and 
the instantaneous rate on beam m, respectively. 

A (r) denotes the average number of users feeding back per beam. It measures the performance of the 
system along the feedback dimension. We have A(t) = Yll=i ^ ^i} since beams are statistically 

identical. We let pi = Pr{7i^i > Tj}. Given a finite feedback system with a feedback constraint A, the 
main optimization problem to be solved is to determine the optimal threshold values to maximize the 
downlink throughput, which is formally written as 

maximize R (r) 

_ (2) 

subject to Y.l=i Pr H,i >Ti}<\ 
As shown in ifTOl . this optimization problem is not easy to solve even for a two-user system due 
to the non-convex objective function and the non-convex constraint set depending on the distribution 
of the SINR values. The complexity of the problem increases with increasing numbers of users due 
to the dimensionality growth. In this paper, we will use the theory of majorization |[T2l to analyze the 
underlying Schur-concave structure in the objective rate function, and solve this optimization problem 
for a general n-user system. Schur-convex/concave functions frequently arise in mathematical analysis 
and engineering applications |[l3l. For example, every function that is convex and symmetric is also a 
Schur-convex function. In the next section, we will briefly introduce some key concepts from the theory 
of majorization to be used in our analysis later. 

III. Majorization 

For a vector x in M", we denote its ordered coordinates by > ■ ■ ■ > For x and y in M", we 
say X majorizes y and write it as x y if we have X]i=i ^(j) — Y^l=i V(i) when j = 1, • • • , n — 1, 
and Y^^=i^{i) — Y^l=iy{i)- ^ function (/? : R" i— M is said to be Schur-convex if x y implies 
ip{x) > ip{y), and if is Schur-concave if —ip is Schur-convex. A Schur-concave function tends to 
increase when the components of its argument become more similar. We will establish conditions under 
which R (t) becomes a Schur-concave function, which will, in turn, imply the optimality of homogenous 
threshold feedback policies. 

However, establishing these conditions is not straightforward, which makes the following lemma vital 
for our analysis. 

Lemma 1: Let ip he a real-valued function defined on M", and D = {z € ffi" '■ zi > Z2 > ■ ■ ■ > Zn}. 
Then, (/? is a Schur-convex function if and only if, for all z € I? and i = 1,2, ■ ■ ■ , n—1, ip{zi, . . . , Zi^i, Zi + e, Zj+i — e, ZiJ^ 
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is increasing in e over the region < e < min (zj-i — Zi, Zi+i — Zj+2)- ' 

It can be seen that the coordinates Zi and Zj+i are systematically altered by using the parameter e, and 
the constraints on e eliminate any violation in the order. Interested readers are referred to [12] for more 
insights on the theory of majorization. In the next section, we will see how we can use this theory to 
identify the Schur-concave structure in the objective rate function. 

IV. SCHUR-CONCAVITY ANALYSIS 

In this section, we will establish sufficient conditions on the SINR distributions for the Schur-concavity 
of the rate function. We focus on the first beam to explain our proof ideas without any loss of generality 
since all beams are statistically identical. We start our analysis by analyzing the rate as a function of 
thresholds and establishing three important lemmas. Then, we will incorporate the feedback constraint 
into our optimization problem by interpreting the rate as a function of feedback probabilities. Using these 
results, we will finally establish the underlying Schur-concave structure in the rate function through the 
theory of majorization. 

A. Rate as a Function of Thresholds 

Consider thresholds in increasing order, i.e., t^[i) < • • • < T^(i) < T^[i+i) • • • < t-k^u)- Now, by analyz- 
ing (r^(i+i) + e, r^(i) - e) = (r^(„), • • • , r^(i+i) + e, r^(i) - e, ■ ■ ■ , r^(i)), we can use Lemma [U 
to identify the underlying Schur-concave structure in the rate function. However, analysis of this function 
is still complex. Therefore, we use the following method to reduce problem complexity. 

Let A/"' = {A: : k e M k k ^ 'i:(i),T:(i + 1)}. We fix the thresholds and the SINR values of all users 
in M' . Randomness is now associated only with users 7r(i) and 7r(i + 1). Let 7_^, = max7jt^il|^^^ i>Tk}- 
The instantaneous rate on beam 1 as a function of 7,r(j).i' 77r(i+i),i ™d 7_^, is (77r(i),ij 77r(i+i),i) Ta/"') ~ 
log (l + max{%(i)_i,%(i+i)_i,7_;^,}). Therefore, 

It is not hard to see that (''"7r(i)7 ''"7r(i+i) ITaA') depends on the value of 7j^,. jp, can be greater than 
Tn(i+i)^ less than T^(j) or in between r^(j) and r^i^i^i). We will now establish three important lemmas for 
these three cases, which will be useful in interpreting the rate. We will start with the case jp, > r^(j+i). 

'At the end points i = 1, i — n ~ I, the condition is modified accordingly. 
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Lemma 2: If 7^, > T^(i+i), (r^(i), r^(i+i) ItJv') only depends on 7^,, and given by G(7j^,). where 
G(7XrO = Pr {7- .+1 < 7Xr'l7>} log (1 + 7>') 

+ E[log(l+7.W)l{7.W: 

and 7*i+i = max {7^(i), 7^(^+1) }. 

Proo/; Let 7*.^! = max {7^(j), 7^(^+1)}. From ©, 

(''"7r{j),T7r{j+l)l7jV') 

= log(l + 7»Pr{7-m<7Xr'l7>} 



log(l + 7-m)l{7.w>7;^„}l^^' 



• (4) 



+ E 

For the first term on the right hand side of ©, it is not hard to show that |7j*j+i < 7j\/-'| — 

< 7jV'} since j^, > r^^i+iy 
For the second term, we must have 7*^4.1 = 7*j+i since 7*j^_i > 7^^-, > T^(i+i), ■ Thus, we end up 
with (r^(i), r^(i+i) |7>') = G{^Xf')- ■ 
The next lemma will evaluate the value of R^ (T7r(i), |7_/^/) for j^, < T^i^iy 

Lemma 3: If j^, < T^^i), R^ {T^(i),T^{i+i)\7^,) is given by 



Ri{Tn{i),'rn(i+i)\lAf') = I log{l+x)dF (x) 

+ F(r^(i+i)) /" 'log(l+x)dF(x) 



Proof: When 7;^, < r^(i), © simplifies to R^ (T^(o,T-^(i+i)l7jv') = 

E [log (1 + 7*i+i) l{7^(,,,i>r,(,),7,(,+i).i>r^(.+i)}] 

+ E [log (1 + 7^{i+i),i) l{7,(.+i,,i>-.(.+i)}] ^ Kw) 
+ E [log (1 + 7,(i),i) 1{^^,.,,,>.,(„}] F (T,(i+i)) 

+ log (1 + 7>,) ^ (^^(i)) ^ (-^-(i+i)) • 

The first three terms on the right hand side is identical to the rate expression of the two-user system 
analyzed in lITOl . Substituting the result for the two-user case completes the proof. ■ 
Finally, we look at the case where r7r(j) < 7^^, < r^(i+i)- 
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P-K{i+l) - mm -p,r(i+2),P7r(j-l) -P7r(i)j 



Fig. 1. Ordered feedback probabilities, and the range of q and e. 



Lemma 4: If T^(j) < 7^, < r^(i+i), then the value of {T^ii),T^(i+i)\llf) is given by 

/•oo 

R2 {r-n{i).'T-K{i+i)\lN') = I log(l + x)dF'^{x) 




log(l + x)dF{x) 



+ log (1 + 7j^,) F (t,(,+i)) (1 - Pr {7,(,) > 7>|7>}) ■ 

This lemma can be proved by using similar arguments to the ones used in the previous two lemmas. 
We skip its proof to avoid repetitions. 

If 7j^, = r^(j), Ri and R2 in Lemmas |3] and |4] evaluate to the same expression. Similarly, if 7j^, = 
T^(i+i), G and R2 in Lemmas |2] and |4] evaluate to the same expression. This shows that the rate as a 
function of 7j^, is continuous at T^(j) and T^(^iJ^i). 

Given the initial threshold values {T^(j)}J_^, the first approach to discover the Schur-concave structure 
in the rate function is to analyze the behavior of the function 

grie) = R^ (T^(i) - e,T^(j+i) + e\^'l^,) 

for e € [0, mill (r^(j) — 'T-K{i-i)^'T-K{i+2) ~ ''"7r(i+i))] by making use of Lemma [T] This is now a scalar 
problem. However, we find it more useful to interpret the rate as a function of feedback probabilities 
since the feedback constraint in Q is in terms of these probabilities. Through this interpretation, we can 
incorporate the feedback constraint into our optimization problem more easily, as will be shown in the 
next sub-section. 

B. Rate as a Function of Feedback Probabilities 

Since r^(j) = F~^{1 — ^^^(i)) and / has the support M+, we can write R^^ {T^{i),T^[i+i)\l}j'') as 
R^ {P-K(:i)TP-K{i+i)\l^f') without any ambiguity. Since F is monotone increasing, we have P7r(i) ^ P-k(2) ^ 
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■ ■ ■ > P-K(i) > P-K{i+i) > ■ ■ ■ > P-K{n)- Therefore, given the initial feedback probabiUties {p7r(j)}J^i> we 
analyze the behavior of the function 

5p(e) = -R^ {P-K{i) + e,P7r(i+i) - e!7» (5) 

for e G [0,min (^^^(j^i) — P-n(i),P-n(i+i) — P7r(j+2))] to discover Schur-concavity of the rate function by 
Lemma [T] and P7r(i+i) give us the feedback level Aj = ^^^(i) + P7r(j+i)' ^rid other probabilities 

give us natural boundaries on ^^^{j) and p-n(i+i) as Pn{i+2) < ?'7r(j+i) < P7r(j) < ?'7r(i-i)- Without violating 
these boundaries, we can vary P7r(i+i) by keeping Aj constant. That is, we can take any q G "Pj+i and write 

(PttW +e,?'7r(i+i) -elTAT') as a function of g, where P^+i = - min - p^(i+2),P7r(i-i) -P7r(i)) , t] ■ 

Please refer to Fig. [T] for a graphical representation of this selection of q. 

Now, by using LemmaHOandU we have i?i {q\TM') = G {T^') ^{g>i-F(7j^,)}+^i (^l7>') '^{q>K-{i-F{ru,))] + 
R2 {q\rM') 1{,<1-F(7;„) & q<x.-{i~F{y;„))} for 9 S Vi+i. 

Some insights about this rate expression are as follows. Assume 1 — F (j^,) < ^. Then, when q 

changes from - min - Pn{i+2),P7r{i~i) - P7r(i)) to the threshold of the user with 

index Tr{i + 1) can become as small as F^^ (l — and can become as large as T^(j_|_2)- Therefore, 
when T^(j_|_i) is large, 7j^, lies in between two thresholds, and the rate expression is given by i?2 (^ItJV') 
for small values of q. When r^(i+2) is small, 7j^, always lies to the right of T^(^i+i), and the rate is given 
by the constant G (7^/). Similar explanations can be given for the rate expression for 1 — F (7j^/) > ^■ 

C. Schur-concavity of the Rate Function 

Now, we show how we can utilize the theory of majorization to identify the Schur-concavity of the 
rate function. First, we show that (q'It'Jv'') an increasing function of q over Vi+i. Note that we are 
proving that the rate function is increasing over a larger region than the one needed in Lemma [T] (Please 
refer to Fig. [T]) The following lemma formally states sufficient conditions for this property to be true. 

Lemma 5: R^ {lll^f'^ increasing function of q for any given jp, if / is bounded at zero, and 
has the derivative /' satisfying /' (^F^^{q)^ < — ({^-i^^j) for all q G 'Pj+i and i = 1, . . . , n — 1. 
Proof: We can write Ri {q\^^,) and R2 {q\^^,) explicitly as 

POD 

RiilllM') = / log{l + x)dF\x) 

JF-^{l~q) 

rF-^il-q) 

+ il-q) / log{l + x)dF{x) 

JF-^l+q-Xi) 

+ log(l + 7j^,)(l-^?)^(l + ^/-Ai), 
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and 



R2{q\lM')= / log(l + x)dF2(x 

JF-^{l-q) 

fF-'{l-q) 

+ {l-q) / log(l + x)dF{x 



+ log(l + 7>0(l-'Z)i^(7>'), 

respectively. We can show that ^^ite) > for any q < ^ by directly using the two-user analysis 
in ifTOl . given the conditions on the distribution are fulfilled. Integration by parts on the derivative of 
(g|7>,) gives us = J^^d"^) ^dx > 0. Since G(7>,) is independent of q, '.EiMA > 

over "Pj+i, which completes the proof. ■ 

The above lemma proves that if the conditions on the SINR distribution are satisfied, {<l\l%-i) is 
an increasing function of q. This implies that R^ {Pn{i) + ^jPn{i+i) ~ ^\lXf') is ^ decreasing function 
over the region of interest for e. (Please refer to Fig. [T]) Therefore, we can claim that the rate is a 
Schur-concave function of the feedback probabilities if the conditions in Lemma [5] are satisfied. This is 
formally stated and proven in the next theorem. 

Theorem 1: R{p) is Schur-concave if / is bounded at zero, and has the derivative /' satisfying 
/' {F-\x)) < - (ftS^ for all x G [0, 1]. 

Proof: By Lemma [51 —R^ {Pn{i) + ^,Pn{i+i) — ^ITaT') increasing function for any given 7j^, 
over e G [O, min — Pn{i),Pn{i+i) ~P7r(i+2))]- Therefore, by using Lemma [T] and taking the 

expectation over 7j^,, we have 

-R^ {Pwii) + e,P^(^i + l) - e) > -R^ {P7T{i),P7T{i + l)) 

for all e G [0,min - P-jrii) , P-Kii+i) -P7r(i+2))] - Since this is correct for all p G [0, 1]", -R^ (p) 

is a Schur-convex function of p, implying that R^ (p) is Schur-concave. Since beams are identical, we 
also conclude that R (p) is a Schur-concave function of feedback probabilities if the conditions in the 
theorem are met. ■ 
A Schur-concave function increases when the dispersion among the components of its argument 
decreases. Therefore, this theorem implies that a solution for the optimization problem in dD is r* = 
{F--^l-^),F-^l-j,),...,F~^{l-j,)yiff'{F-Hx)) < -((g^ for all X G [0,1]. 
Next, we will see how our results can be applied to a Rayleigh fading channel model, which is one of 
the well known channel models in the literature, e.g., ||T41 . |[T5l . and closely approximates measured data 
rates in densely populated urban areas [|16l . 
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V. Applications: Rayleigh Fading Channels 

In this example, we will assume that transmitted signals and Rayleigh distributed channel fading coeffi 
cients are of unit power. Interested readers are referred to |3J for further details of the physical layer model 
In this case, F and / can be given as F{x) = 1 — (^^^-^yi-i and f{x) = 



+ + M- 1 



and the functional inverse of F, F ^, is given as F ^{x) = — plog (1 — x) if M = 1 and F ^(x) 



-1 + (M - l)pW - x)^ J if M > 2, where x G [0, 1] and W is the Lambert W 

function given by the defining equation W{x) exjp{W{x)) = x for x > —K 
We summarize our findings for Rayleigh fading environments as follows. 

Theorem 2: For a Rayleigh fading environment with M = 1, R {p) is a Schur-concave function of 
feedback probabilities, and therefore a homogenous threshold feedback policy solves Q if p < 1. 

According to this theorem, a homogenous threshold feedback policy will be optimal if the background 
noise power is large enough, i.e., signal-to-noise-ratio less than OdB. Although, intuitively, a homogenous 
threshold feedback policy may sound optimal for all channel conditions, a direct counter example can 
be given for p > 1 and M = 1 [10]. 

However, our next theorem is promising for Rayleigh fading channels, and it shows that, homogenous 
feedback policies are optimal for a broad range of model parameters. 

Theorem 3: For a Rayleigh fading environment, R (p) is a Schur-concave function of feedback prob- 
abilities, and therefore a homogenous threshold feedback policy solves Q if M > 1. 

Both theorems are direct applications of Theorem [T] Since the multi-beam scenario is a more general 
case, one may a priori assume that homogenous feedback policies should always be optimal for M = 1 
since they are optimal for M > 1. However, this intuition is not correct, and homogenous threshold 
optimality for the single beam case depends on the background noise level. Some more practical insights 
on this result can be found in |10|. 

VI. Conclusions 

In this paper, we have analyzed optimal threshold selection problem for threshold feedback policies 
for finite feedback systems. This is a non-convex optimization problem. We have established sufficient 
conditions under which the rate function becomes a Schur-concave function by utilizing the theory of 
majorization. Schur-concavity of the rate function automatically implies the optimality of homogenous 
threshold feedback policies in which all users use the same threshold value for their feedback decisions. 
We illustrate an example application of our results to Rayleigh fading environments. More refined 
sufficient conditions on the number of beams and the background noise level for the homogenous threshold 
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optimality have been obtained in this case. It is also seen that homogenous threshold feedback pohcies 
may not be optimal for a single beam system with low background noise levels. 
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